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ABSTRACT 

Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic 
and dynamical properties of early-type galaxies. However, no observations have yet been able to 
disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the 
CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. 

We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy 
towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, 
naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If 
we assume that - within the star formation environment of early-type galaxies - metallicity is the main 
driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is 
necessary. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: fundamental parameters - 
galaxies: stellar content - galaxies: elliptical 


1. INTRODUCTION 

The stellar initial mass function (IMF) determines the 
ratio between dwarf and giant stars, and therefore, ulti¬ 
mately regulates the chemical enrichment and the stel¬ 
lar feedback. Moreover, it is the Rosetta stone in our 
understanding of unresolved stellar systems, relating the 
properties of resolved nearby stars to the integrated light 
of more distant objects. 

Over the last years, there is growing evidence of a vari¬ 
able IMF, as opposed to the common as sumption that 
the IMF of the Mi lky Way is universal (|Kroupal l2002t 
iBastian et all 12010 1. These claims come from a wide 
varie ty of approaches, inc l uding stellar population anal¬ 
ysis dCenarro et~ahl 120031 Ivan Dokkurn fa Comm! [20 1 Ol : 
iSpiniello et al. I2012t iFerreras et al.l 120131), gra vitational 
tensing dTreu et aT~ 1 201 Ofc Thomas et all I201T1) and dy¬ 
namical models ( Cappellari et al.ll2012f) . 

Assuming that the answer to all of these observational 
results is indeed a non-universal IMF, two main param¬ 
eters have been proposed. On the one hand, the inte¬ 
grated stellar velocity dispersion is found to correlate 
with the I MF slope as inferred fr om both dynamical 
models (e.g. ICappcllar i et al.l 120121) and stellar popula¬ 
tion a nalysis (e.g. La Barbera et a l. 2013T) . On the other 
hand, iConrov & van Dokkuml (|2012 ) claimed a stronger 
correlation when comparing the [Mg/Fe] pattern to the 
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IMF slope. Nevertheless, these trends have to be care¬ 
fully interpreted because they are based on the spatially- 
integrated light of early-type galaxies (ETCs). Thus, the 
inferred properties are luminosity weighted quantities, 
convolved by the radial light distribution, which depends 
itself on the stellar velocity dispersion (Graham et al. 
120011) . and is hi ghly peaked towards the centr e. 

We sh o wed in Martfn-Navarro et al. I (f2015alfbl. see also 
iPastorello et al.l 20141) that the IMF is a local property, 
not only varying among but also within galaxies. The 
lack of a large statistical sample, however, limited further 
interpretations of our data, and the question about which 
is the main driver of IMF variations remains open. In this 
Letter, we use the 2D spatially-res olved IFU data from 
the CALIFA survey (jSanchez et ablliioTiil 'l to address this 
problem by analyzing the radial IMF profiles of a sample 
of 24 ETGs. 


2. DATA 

We selected the sub-sample of 24 ETGs with redshift 
0.018 < z < 0.030 among the o bserved C ALIFA galaxies 
(for sample properties see iWalcher et al.l 120141) . In the 
local Universe, the strong IMF sensitive index Ti02 is 
potentially affected by a telluric absorption feature at 
~6280A. The selected redshift window maximizes the 
number of ETGs unaffected by this feat ure. Technical 
detail s of the data are fully descri bed in iSanchez et al.1 
(|2012f) : iGarcia-Benito et all (1 201 5fl . The central veloc- 
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Fig. 1.— The best-fitting IMF slope is compared to the local a (a), Vi- ms (b), [Mg/Fe] (c) and age (b). Neither the kinematics 
properties nor the [Mg/Fe] or the age follow the measured IMF variations (p a = 0.35, py rms = 0.30, P[Mg/Fe] = 0.21, p age = —0.50, with p 
being the Spearman correlation coefficient). The right vertical axis represents the IMF slope in terms of F 0 . 5 , defined as the fraction (with 
respect to the total mass) of stars with masses below 0.5 M©. 


ity dispersion in our galaxy sample ranges from ~160 
to ~310 kins -1 , although the bulk of them (21 galax¬ 
ies) have velocity dispersions greater than 200 kms -1 . 
The mean st ellar mass in our sam ple is A£* = 10 1154 Mq 
according to lGonzalez Delgado et al.l (I2014H . 

We derived the line-of-sight mean stellar velocity (V) 
and velocity dispersion (a) following Falcon-Barroso 
et al. (2015, in prep.) Using these measurements, 
each IFU spaxel was then corrected to the rest-frame. 
Finally, we radially binned each galaxy using elliptical 
apertures. The size of the apertures was set to reach 
a signal-to-noise of, at least, 125 per A at AA = 6000, 
6200A. All the quantities discussed in this Letter are 
averaged spaxel-values over these elliptical annuli. 


3. ANALYSIS 
3.1. Stellar populations 

For this work, we made use of the exten ded M I LES 
(MIU SCAT) stellar population models (lYazdekis et al.1 
mi2tm We assume a bimodal, low-mass tapered, 
IMF, whose only free parameter, Tb, is the slope of 
the high-mass end (above 0.6M Q ) of t he distri b ution. 
This p arametrization, first introduced bv lVazdekis et ahl 
ljl996l) . generalizes the Kroupa IMF, which is recovered 
for Tb = 1.35. Note that the current version of the ex¬ 
tended MILES models cover from Tb = 0.3 to Tb = 3.3, 
which can lead to some saturation for high IMF values. 

Given the CALIFA wavelength range 
(AA 3700,7500 A), we focused on five pr ominent spectr al 
indic es: H^ 0 , [MgFe]', Mg2Fe dBruzual fc Chariot! 

I2003T) . NaD, Ti02 CALIFA and TiOi. The last three are 
IMF-sensitive features, whereas [MgFe]' and Mg2Fe 
depend on metallicity and H^ 0 on age. Th e TiO a^ 
index follows the standard Ti02 definition (iTraeer et al.1 
I1998H but with a narrower blue pseudo-continuum 
(AA 6060.625, 6080.625 A) to avoid any telluric contam¬ 
ination. The correction of H A, from nebular emissio n 
was done in the same way as in lLa Barbera et al.1 (|2013t l. 
For each spectrum (i.e. each galaxy and radial bin), we 
inferred stellar population properties with three different 
approaches: 


i) Following lLa Barbera et all (120131) . we minimized: 


X 2 (F b , age, [M/H]) 



A a ,i) — EWm,i 

&EWi 


(!) 

where our index measurements (EW t ). after being cor¬ 
rected for non-solar abundances (A Q) j), are compared to 
the model predictions ( EWm,i )■ The fitting process was 
repeated by considering different combinations of IMF 
sensitive indices, leading to consistent results. 

ii) We implemented an iterative fitting scheme where 
age and metallicity were calculated first using an index- 
index (Hg 0 -[MgFe]') diagram, assuming a standard IMF. 
The age estimate coming from this first step was then 
used to derive the metallicity and the IMF in a second 
index-index (Ti02 CALIPA [MgFe] / ) diagram. These two 
steps were repeated, re-deriving age and metallicity from 
the Hg 0 -[MgFe]' diagram with an updated IMF, until the 
solution converged. 

iii) Alternatively, we also fit our data following Eq. [T] 
but neglecting individual [a/Fe] corrections. 

The three approaches show a good agreement, with 
small differences only for low IMF-slope values, where 
IMF effects can be mimicked by changes in other stel¬ 
lar population parameters. For simplicity, we will re¬ 
fer through this Letter to the best-fitting values derived 
from the simultaneous fit (Eq. [L} of the H^q, [MgFe]' 
and Ti02 CALIPA indices. Note that the latter, com¬ 
pared to the TiOi, exhibits a milder depende nce o n both 
abun dance ratios and total metallicity (lLa Barbera et ahl 
l2013f) . However, the absolute Ti02 sensitivity to these 
parameters dep ends on the adopted stellar population 
synth esis model ([Bruzual fc Gharlotll2003i : iTlromas et all 
2011). Moreover, Ti02 is less sensitiv e to variations 
in the temperature scale of giant stars. Spi niello et al.1 
(|2015f) showed that, after accounting for metallicity, no 
variation in the effective temperature is needed to fit the 
strength of gravity-sensitive features in massive ETGs, 
and that a similar result would be recovered with other 
stellar population models if the same IMF parameteriza¬ 
tion were used. Although our fitting scheme does not ac¬ 
count for the residual impact of non-solar abundances on 
the temperature scale of the (solar-scaled) isochrones, the 
[a/Fe] effect on the isochrones is si gnificantly milder tha n 
that on the stellar atmospheres llVazdekis et al.l [2015T) . 
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Fig. 2.— IMF-metallic ity relation obtained f r om CA LIFA local measurements (blue). We also show the local IMF and metallicities 
measurements derived by I Mart in-Navarro et ahl J2015al ibT) (red, orange) for three of nearby ETGs, as well as global SDSS measurements 
(black). We found it to be the strongest correlation (p[m/h] — 0.82). As in Fig. [T] the right vertical axis indicates the F 0.5 ratio. For 
reference, the standard Kroupa IMF value is shown as a horizontal dotted line. Dashed line correspond to the best-fitting linear relation 
to all the datasets. 


The latter is corrected in our approach (A a * in Eq. [TJ . 

In addition to the IMF, age and metalliciy, the [Mg/Fe] 
of each radial bin was derived b y means of the [ZMg/Zpe] 
proxy (lLa Barbera et al .11 20131 ). i.e., using the metallic- 
ity difference between two index-index diagrams, where 
H/ 3 0 is plotted against a Mg and Fe metallicity indicator, 
respectively. 

3.2. Stellar kinematics 

To understand the IMF variations, we also compared 
our best-fitting IMF values to two kinematics parame¬ 
ters: the local a and the local H rms defined as H rms = 
y/V 2 + tr 2 . Whereas a has been clai med t o be the 
main driver of the IMF variations dTreu et all 120101 
iFerreras et al.l I2013T ) . in spatially resolved studies V-ms 
accounts for both random and ordered motions. 

4. RESULTS 

In Fig. |T] we present the correlation between the best¬ 
fitting IMF slope, Tb, and the local values of [Mg/Fe], 
age, <J and I/. ms . None of them show a tight correlation 
with Tb- The mild relation between age and r b can be 
understood either as a residual degeneracy between both 
parameters or as a consequence of the IMF-metallicity re¬ 
lation, since young stars within massive ETGs are likely 
formed from metal-enriched material. 

4.1. The IMF-metallicity relation 

Among all the explored relations, the IMF slope- 
local metallicity relation emerges as the most funda¬ 
mental. This is shown in Fig. [2j where the local r b - 
[M/H] relation derived from the CALIFA survey (blue) 
is combined with local IMF estimates, obtained at dif- 
fer ent ga l actocentric distances, for three nearby ETGs 
bv lMartin-Navarro et all (l2015al lbl) (red and orange sym¬ 
bols). In addition, we also show (black) the best-fitting 
IMF and metallicity inferred from SDS S stacked spectra. 
To con struct these spectra, we followed lLa Barbera et all 
(12013T) . but binning according to both a and [M/H] of 
the individual galaxies. The broad wavelength range 


of the SDSS data set allows us to infer the IMF, not 
only using those features within the CALIFA spectral 
range, but also prominent near-IR IMF sensitive fea- 
tu res such as the NaI8189 and the Call triplet (see §4.1 
in lLa Barbera et al.l 120131 for a detailed description of 
gravity-sensitive features in SDSS spectra). 

The fact that the three datasets included, although 
based on different sets of line-strengths, lie on the same 
relation, supports a tight connection between IMF slope 
and metallicity, regardless of the details in the determi¬ 
nation of the stellar population parameters. Moreover, 
the agreement between integrated measurements from 
the SDSS spectra and spatially-resolved values, suggests 
that the mechanism behind the local IMF variations ul¬ 
timately shapes the global galaxy mass-IMF relation. 

A linear fit to all the measurements shown in Fig. [2] 
leads to the following relation between IMF slope and 
metallicity in ETGs 

r b = 2.2(±0.1) + 3.1(±0.5) x [M/H] (2) 

Since IMF-sensitive features ul timately trace the dwarf - 
to-giant ratio F 0 . 5 , as defined in lLa Barbera et all (1201311 . 
the above equation can be expressed in terms of a single 
power law IMF as 

T = 1.50(±0.05) + 2.1(±0.2) x [M/H] (3) 

Apart from the measurement errors, the scatter in the 
relation comes from two sources: the IMF-[a/Fe] de¬ 
generacy when fitting gravity-s ensitive features around 
the Kroupa-like IMF regime (lLa Barbera et ahl 12013T) 
and the dep endence o f the I MF on the minimized set 
of indicetQ (|Spiniello et al.l l2014af ) . In this sense, the 
fact that the Ti02-based CALIFA measurements show 
a steeper IMF-metallicity trend is consistent with a 
stronger metallicity dependence of this index than pre¬ 
dicted in the MILES models. On the other hand, the con¬ 
sistency among different data-sets (p = 0.82 when CAL- 

1 Uncertainties in Eq. 2-3 account for this effect, by repeating 
the fit using only CALIFA data. 
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[MgFe]' 


Fig. 3. — Empirical relation between the metallicity-sensitive 
[MgEe]' and the IMF-sensitive Ti02 CALIFA features. Index mea¬ 
surements (at a resolution of 200 kms -1 ) are color-coded by their 
value, as an age proxy. An IMF-metallicity relation is needed 
to explain the observed trend, since the Ti02 CALIFA weakly de¬ 
pends on the total metallicity and [MgFe]' is almost independent 
of the IMF. 

IFA, SDSS and IMartin-Navarro et al.l (l2015al lbl) are con¬ 
sidered) points to a genuine IMF-[M/H] trend, as shown 
in Fig. [2] 

4.2. The [MgFef-Ti02 CAljlFA empirical relation 

To strengthen the validity of our result, we adopted an 
empirical approach. In Fig. [3] we compare the [MgFe]' 
to the Ti02 CALIFA individual measurements for all radial 
bins in our sample, at a common 200 kms -1 resolution. 
Each point is color-coded by its H^ 0 value, represent¬ 
ing an age segregation. The [MgFe]' index is indepen¬ 
dent of the IMF, depending almost entirely on the total 
metallicity. On the contrary, the TiOar,., ,^ is a strong 
IMF indicator, that increases with age (lLa Barbera et all 
iM3). Since MILES stellar population models predict 
the Ti02 CALIPA feature to be, in the explored metallicity 
and age regime, almost independent of total metallic- 
ity, the strong (correlation coefficient p = 0.86) correla¬ 
tion shown in Fig. [3] can be interpreted as a metallicity- 
IMF t rend. A similar relation was found bv lTraeer et al.l 
(11998T ) while comparing the Mg 2 and the Ti02 indices. 
Note that if the Ti02 sensitivity to total metallicity is 
larger than predicted by MILES models, a certain corre¬ 
lation between these two indices would also be expected. 

5. DISCUSSION 

ETGs are characterized by old stellar populations and 
thus, only stars with masses M < 1M 0 remain alive. 
Therefore, the analysis of gravity-sensitive features in 
the integrated light of unresolved ETGs can only con¬ 
straint the the dwarf-to-giant ratio. We vary this ratio 
by changing the slope (Tb) of the high-mass end of the 
IMF, while it was kept constant for stars with masses 
below 0.5M e . Notice that our results, restricted to in¬ 
ference of the dwarf-to-giant ratio, barely depend on the 
detailed IMF shape, but it must be considered when ex¬ 
ploring, for example, the chemical evolution and the ex¬ 
pected mass-to-light ratios of ETGs. 

5.1. The underlying parameters behind the varying 
dwarf-to-giant ratio in ETGs 


Two competing candidates have been proposed to 
explain the observed dwarf- to-gian t ratio variations in 
ETGs: <7 and [Mg/Fe]. ISmithl (1201411 investigated 
these two parameters o v er the same sample of galaxie s 
(iGannellari et al.l 12011 : iGonrov fo van Dokkumi I2012H . 
pointing out that the stellar population analysis favors 
the stellar population property ([Mg/Fe]), whereas a dy¬ 
namical analysis supports the dynamical-related quan¬ 
tity (cr). The work of iSmithl opened the question of 
whether IMF studies are actually probing the IMF or 
whether they are, at least partially, dominated by con¬ 
founding factors. However, s tellar population analy¬ 
ses of SDSS stacked spectra (|La Barbera et al.l 120131 : 
iSpiniello et al.l l2014af) have shown a strong correlation 
between a and t he dw arf-to-giant ratio. Moreover, 
lLa Barbera et al.l (|2015H . using also SDSS stacked spec¬ 
tra, showed that [Mg/Fe] is loosely (or un-) correlated 
with the IMF slope. The latter result is in agreement 
with panel (c) in Fig. [I] where we show that the local 
[Mg/Fe] is decoupled from the dwarf-to-giant ratio vari¬ 
ations. 

Regarding a, in iMartin-Navarro et al.l ( 2015b ) we 
showed that its local value was not the main driver be¬ 
hind the dwarf-to-giant ratio variations of the massive 
relic galaxy NGC 1277. In the present work we con¬ 
firm this result using a statistically significant sample 
of ETGs (panel (a) in Fig. [I]). Notice that previous 
studies suggesting a connection between a and the IMF 
slope were based on spatially unresolved spectra. Thus, 
their velocity dispersions trace the overall galaxy poten¬ 
tial (mass) rather than the detailed kinematics. Thanks 
to the CALIFA data, we propose that this connection 
between galaxy mass and IMF slope arises from consid¬ 
ering simultaneously both the galaxy mass-metallicity 
and the metallicity-IMF slope relations. In this sense, 
t he dwarf-to- gia nt rat io gradients observed in ETGs 
(IMartin-Navarro et al.1 l2015al lbl) can be also accounted 
for the radial metallicity variation within these galax¬ 
ies. To illustrate the suggested metallicity-driven galaxy 
mass-IMF slope relation, we used the above described 
Sloan data (black dots in Fig. 0 to calculate the mean 
global metallicity per galaxy mass (cr) bin. Then, we 
transformed this metallicity into the expected IMF slope 
according to Eq. 0 The result is summarized in Fig. 0 
Although it can not be directly compared with previous 
works, we qualitatively recover the observed relation be¬ 
tween galaxy mass and IMF slope, with more massive 
galaxies being bottom-heavier. 

5.2. Metallicity as a driver of IMF variations 

Metallicity-driven dwarf-to-giant ratio variations can 
be explained if molecular clouds collapse following a 
Jeans spectrum. Such a scenario is expected if ther¬ 
mal pressure dominates the process, although a turbu¬ 
lent medium can y ield similar results if dust cooling is 
significant dLarsonl 1200511 . Perhaps, a realistic scenario 
would couple a Jeans-driven mechani sm with turbulence- 
induced fragmentation (e.g. IPadoan fc Nordlundl [2002 1. 
Nevertheless, this letter shows the importance that 
metallicity plays in determining the stellar mass spec¬ 
trum. Such a scenario is consistent with studies 
of nearby resol ved gala xies dGeha et al.l I2013D . line- 
strength indices (fCenarro et al.ll2003D and optical colors 
(|Ricciardelli et all 20121) . In low-metallicitv regions, it is 
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found that molecular clouds are very efficient in form¬ 
ing massive stars, leading to top-heavy IMFs. This is 
supported by observations of metal-poor systems show¬ 
ing ev i dence of a n IMF dominated by massive stars 
(IMarks et all 120 lH ). 

If metallicity is one of the main drivers behind the IMF 
variations, a number of fundamental aspects of galaxy 
evolution must be revisited. Because the stellar feedback 
regulates the chemical enrichment of galaxies, IMF, star 
formation and chemical composition will be tightly re¬ 
lated during the time evolution. In particular, massive 
galaxies are expected to increase their metallicity dur¬ 
ing the formation of their stellar populations, and there¬ 
fore, the dwarf-to-giant ratio should have been smaller 
at higher redshifts, such that there were fewer low- 
mass stars per high-mass star. Such a time-dependent 


IMF scenario (lArnaud et al 


ILarsonl Il998t IWeidner et all 


1992; 


2013c 


Vazdekis et al 


Ferrer as et al. 


■ 

1 120151) 


would explain the observations of nearby massive ETGs, 
which are inconsistent with a ti me-invariant, steep ( or 
even Milky-Way like) IMF slope (lArrigoni et alJl2010D . 

Moreover, the [Mg/Fe] has been extens ively used as a 
proxy for the formation time-scale (e.g. iThomas et all 
12005 ?): rapid formation events lead to high [Mg/Fe] 
values. However, an IMF dominated at early times 
by high-mass stars would also produce an enhanced 
[Mg/Fe], relaxing the typical (~ lGyr) constraint on 
the rapid formation time-scale. Thus, to safely inter¬ 
pret z ~ 0 observations, it is crucial to understand the 
chemical evolution of galaxies since they were formed, 


properties (IVazdekis et al. 1996 

). In this sense, it has 

been shown (e.g. 

Ferreras et al. 

2009: Choi et al. 2014 

Ferrc-Mateu et al. 

2014|) that the chemical composition 


of massive galaxies has remained constant over the last 
~ 7 Gyr. If metallicity actually regulates the dwarf- 
to-giant ratio variations, it would imply that the IMF 
of massive objects was bottom-he avy at z ~ 1, as re¬ 
cently suggested by o bservations dShettv fe Capnellaril 
12014 iMartm-Navarro et all I2015d ). since no chemical 
evolution has happened since then. 

Finally, we want to emphasize that the present study 
is restricted to the analysis of ETGs. The local condi¬ 
tions (e.g. turbulence, pressure, density, radiation, and 
magnetic fields) of the interstellar medium at high- 2 , or 
even in nearby disk galaxies, are expected to differ from 
the local environment within ETGs. Thus, more observa¬ 
tional efforts are still needed to fully characterize the IMF 
behavior with varying star-f orming co ndi tions. In fact, 
for very metal-poor systems, iGeha et all (|2013l) found a 
slightly top-heavy IMF, but steeper than expected from a 
crude extrapolation of Eq.[2] The trend in Fig.[2]seems to 
actually flatten at lowest and highest metallicity, point¬ 
ing to a non-linearity of the IMF-metallicity relation. 


6. SUMMARY 

Although deviations from the standard Milky-Way 
IMF have been extensively reported in the literature 
over the last years, the mechanism responsible for these 


variations remains unknown. We have analysed the ra¬ 
dial gradients of IMF-sensitive features in a sample of 
24 ETGs observed by the CALIFA survey, finding that 
the local IMF is tightly related to the local metallicity. 
Our result agrees with previous works reporting local 
Galaxy mass [log MJ 

10.50 11.00 11.50 12.00 



Galaxy velocity dispersion [km/s] 


Fig. 4. — Predicted IMF slope — galaxy velocity dispersion re¬ 
lation. The best-fitting metallicity of the global SDSS measure¬ 
ments (right vertical axis) was transformed, following Eq. [2] into 
an expected IMF-slope value Fy. The upper horizontal axis rep- 
resents the galaxy ma ss, estimated from the velocity dispersion 
(ICappellari et all 12013 ). We suggest that the relation found be¬ 
tween IMF slope and galaxy mass can be understood as the com¬ 
bination of the mass - metallicity and the metallicity - IMF slope 
relations. 


and global IMF variations, and it explains the observed 
galaxy mass-IMF relation. We speculate about the im¬ 
plications of metallicity-driven IMF variations in the con¬ 
text of galaxy formation and evolution. 

The intimate connection between IMF and metallicity 
described in our work suggests a complex massive galaxy 
formation process, departing from the classical picture 
where these objects formed under nearly stationary 
conditions. Thus, to safely understand and interpret 
2^0 observations, it is necessary to untangle the early 
life of massive galaxies at higli- 2 . The combination of 
more sophisticated stell a r popu l ation synthesis models 
rtGonrov fo van Dokkuml l2012t IVazdekis et, all 120151) 
and high -2 spectr oscopic surveys ( Brammer et al.ll2012t 
iKriek et al .1120141) will provide valuable insights. 
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